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Abstract 

The top quark flavor changing neutral current (FCNC) process is an excellent probe to search 
for new physics in top sector since the Standard Model expectation is extremely suppressed. We 
explore Higgs-mediated top quark FCNC, focusing on H-t-c Yukawa coupling Xct within the general 
two Higgs doublet model. After electroweak symmetry breaking the top quark FCNC couplings 
are included in the charged Higgs Yukawa sector so that they contribute to various processes 
in flavor physics. To probe Act, we study anomalous single top production and the same sign 
top pair production at the LHC in association with flavor physics from the tree-level processes 
B TV, B ^ TV as well as from the loop-level processes B^ —)• B^ g — Bd,s mixing. We 

perform combined analysis of all the constraints regarding the fine-tuning argument to fit the data 
and discuss future prospect. The recently updated measurements on H —)• prefer large 

Xct, but we show that the current bound on the same sign top pair production at the LHC gives 
the most significant upper bound on Xct to be less than 10 ~ 30 depending on neutral heavy Higgs 
masses. We also find that for the given upper bound on Xct, B —)• put significant lower 

bound on H-t-t Yukawa coupling, and the bound is proportional to the charged Higgs mass. 
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I. INTRODUCTION 


The top quark, the heaviest particle in the Standard Model (SM), plays an important 
role as an input for the electroweak (EW) precision measurements [1]. Because its mass is 
much heavier than other known particles, the top quark is considered to be the most viable 
candidate which has a close connection to new physics (NP) that controls the EW symmetry 
breaking mechanism. Meanwhile, the discovery of the SM-like Higgs boson at the LHC [2lE] 
and the precision measurement of its property [H |5] shed much light on the physics in EW 
sector, boosting the relevant studies. Especially, NP scenarios with extended Higgs sector 
have received great interest due to its rich phenomenology and attempt to complement the 
SM [6l [7]. 

One of the simplest scenarios with extended Higgs sector is to introduce a new Higgs dou¬ 
blet. Because the two Higgs doublets can couple to both up-type and down-type quarks, after 
rotating into their mass eigenstates, the tree-level flavor changing neutral current (FCNC) 
inevitably arises. In the SM, the tree-level FCNC is forbidden by the GIM mechanism [8]. 
The FCNC process only takes place through the loop diagrams with charged current and 
rough estimation of the loop correction at the amplitude level is 

where V^km CKM matrices, niq is the mass of quark inside the loop. Thus, the loop- 
induced down-type quark FCNC processes such as 6 —)■ sy, which is involved with top 
quark loop, has enhancement factor {mt/my^Y their rates mostly fall within current 
experimental reach of B physics and Kaon physics. Therefore, the down-type quark FCNC 
is severely constrained and dangerous to many NP scenarios. On the other hand, the up- 
type quark FCNC processes, for example top quark FCNC process t —)■ cy, are involved 
with 6-quark loop and extremely suppressed by {mb/mwY■ The estimation of B(t —)■ cy) is 
[H] within the SM, far too much behind the current experimental reach. 

In order to avoid tree-level FCNC, one usually introduces a discrete Z 2 symmetry to 
make each up-type or down-type quark couple to only one Higgs doublet. In the Minimal 
Supersymmetric Standard Model (MSSM) the supersymmetry itself plays the role. Without 
such a Z 2 symmetry, the general 2HDM which is called “2HDM type HI” follows a specific 
scheme to circumvent severe down-type quark FCNC constraints such as the natural flavor 
conservation HD!, the minimal flavor violation mHis] and Cheng-Sher ansatz HU- In this 
work we adopt the last one, in which the Yukawa coupling connecting quarks with flavor 
indices i,j to one of the neutral Higgses is described as 

= ^ ( 2 ) 

where v is the SM vacuum expectation value (vev), v = 246 GeV, A^- is considered to be 
0{1). With this ansatz, down-type quark FGNG is severely suppressed due to the small 
masses of u, d, s quarks, being safe against the experimental constraints. However, top 
quark FGNG process can be potentially large and should be explored in collider physics as 
well as in flavor physics. 
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In the 2HDM type III, after EW symmetry breaking the top quark FCNC Yukawa cou¬ 
plings Xqt (<? = c) also come into play in charged Higgs Yukawa couplings. Therefore, 
the phenomenology of top quark FCNC process with neutral Higgs exchange is naturally in 
connection with flavor physics process with charged Higgs exchanged due to the common 
Yukawa couplings Xqt- Studies on the top quark FCNC in collider physics especially through 
anomalous top quark decays were performed in Refs. [T8II22] . There have been studies on the 
issue that large top quark FCNC coupling Xct is needed [23l|2l] to explain the measurements 
of B{Bd tv) at BaBar [25], which were quite larger than the SM expectations. The 

authors of Ref. |26| study the collider signature with constraints from 6 —)■ sy concerning the 
perturbativity of Yukawa couplings within the 2HDM and the MSSM. For more comprehen¬ 
sive study on 2HDM type HI contribution to both collider and flavor physics, we refer to 
Ref. [27]. The model independent approach using low energy effective operators was done 
in Ref. [28] . 

In this work we focus on H-t-c FCNC coupling within 2HDM type HI by adopting 
Cheng-Sher ansatz. We perform detailed study on several experimental observables that 
can give bound on X^t from collider physics and flavor physics with the most up-to-date 
experimental data. The issue on B{Bd —)■ D^*'^tv) is revisited with new data from Belle and 
LHCb. Especially it will be shown that the search for the same sign top pair production at 
the LHC plays crucial role to constrain X^- Since the current precision measurements of the 
SM Higgs properties are very well consistent with the SM expectations HE], we assume the 
alignment limit for the Higgs potential of 2HDM type HI, in which the SM Higgs sector is 
well decoupled from the NP sector. 

The paper is organized as follows. In section H, we briefly describe and discuss about 
the Yukawa structure of aligned 2HDM type HI. Section HI explains about the method of 
numerical analysis in this work. In section IV, we study the top quark FCNC processes 
and investigate the bounds from the LHC experiment. In section V and VI, we study the 
constraints from the flavor physics with tree-level and loop-level processes. Section VH 
is reserved for the combined analysis and future prospect for the constraints on X^- We 
conclude and summarize our result in section VIII. 


II. YUKAWA SECTOR OF ALIGNED 2HDM TYPE III 

The Yukawa interaction Lagrangian of 2HDM type HI can be described as [29] 

-Cy = + Ql{Y,^^i + Y^^^ 2 )ur + + h.c., (3) 

where Ql,Lr are left-handed quark and lepton doublets while UR,dR,eji are right-handed 
singlets in interaction basis. The two Higgs doublets <I>i and <I >2 are introduced with the 
dehnition = ia 2 ^* where ct 2 is Pauli matrix. are corresponding Yukawa matrices 

where the flavor indices are implicitly considered. After the EW symmetry breaking $i and 
$2 have the vevs {^i) = Vi/y/2 which satishes vf + V 2 = where v = 246 GeV. As usual, 
we dehne tan/3 = V 2 /V 1 . 
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Then, we diagonalize mass matrices for fermions from Eq. ([^ and for Higgses from 
Higgs potential Lagrangian which is described in many literatures (We refer to review paper 
Ref. Cl)- We define a as a mixing angle of neutral CP-even Higgses. As we discussed in the 
introduction, we adopt the alignment limit that specihes 


sin(/? — a) = 1 , 


(4) 


to make the model comply with the Higgs precision measurement jSDHSZ]- With this align¬ 
ment limit, the Yukawa Lagrangian Eq. (|^ is re-expressed in terms of mass eigenstates as 
follows 


Cy = Cy,sm + ^ d^dH + ^ uCuH + ^ ^ h^^dA - ^ 


+ \u(CVckmPl - VckmCPR] dH+ - uePdH+ + h.c. 


(5) 


by ignoring Goldstone Lagrangian. Here, £y, sm is equal to the SM Yukawa Lagrangian, 
M, d, £ are mass eigenstates of up- and down-type quarks and leptons, H, A are CP-even and 
-odd neutral Higgses, and are charged Higgses. Vckm is the CKM matrix, Pl and Pr 
are chiral projection operators, Pl,r - I (1 = 1 = 75 ). Note that in the alignment limit, the SM 
Yukawa sector is completely decoupled from the NP sector. are Yukawa matrices for 

the mass eigenstates which include all the FCNC couplings. 

In this work we assume that the new Yukawa matrices are CP-conserving, that is 
are real and symmetric : 

^u.d.i /-u.d.i /n\ 

iij = iij = iji ■ ( 6 ) 

To avoid severe constraints from down-type quark FCNC, we adopt Cheng-Sher ansatz, 
Eq (§. Due to the tiny masses of u, d, s quarks, the elements of Yukawa couplings that 
contain those quarks are negligibly small: 


0 0 ^ 


s — 


0 0 e 


sb 


^0 ^sb ^bbJ 






^0 0 
0 
^0 


0 ^ 
^ct 


(7) 


Here, we include ^sb since it can play some role in our study. In this set-up, the only relevant 
top-quark FCNC coupling is Act where ^ct = ■ R should be emphasized that 

the top quark FCNC coupling Act not only belongs to neutral Higgs Yukawa sector but also 
comes into play in charged Higgs Yukawa sector as can be seen in Eq. ([^. This important 
feature leads us to probe Act with the combined analysis of phenomenologies of both collider 
physics via neutral Higgs exchange and flavor physics via charged Higgs exchange. 


III. METHOD OF NUMERICAL ANALYSIS 

Before we discuss the phenomenology of top quark FCNC, we hrst summarize theoretical 
input parameters as well as experimental values that are used in this work and discuss about 
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0.21664 ± 0.00048 

m 

\Vub\ (semi-leptonic) (3.70 ± 0.12 ± 0.26) x 10 ^ 

[38] 

\Vcb\ (semi-leptonic) (41.0 ± 0.33 ± 0.74) x 10“^ 

m 

7r] 

tjo 9+6.3 

0 

[38] 

Wicirnc) 

(1.286 ± 0.013 ± 0.040) GeV [38] 

mb{mb) 

(4.18 ±0.03) GeV 

m 


(165.95 ±0.35 ±0.64) GeV 

[38] 

/f^"(0) 

0.9641 ± 0.0015 ± 0.0045 

m 

JBs 

(225.6 ±1.1 ±5.4) MeV 

[38] 

fBs/fBa 

1.205 ±0.004 ±0.007 

m 

Bb. 

1.320 ±0.017 ±0.030 

[38] 

Bbs/Bb^ 

1.023 ±0.013 ±0.014 

m 


TABLE I. The theoretical input parameters used in the numerical analysis. 


the details of numerical analysis. Table [I] shows input parameters for the processes in flavor 
physics. The values are taken from the latest result of CKMfitter collaboration |39]. To 
obtain the uncertainties of theory prediction, we vary each parameter value within la range 
and add each individual uncertainty in quadrature. 

In Table we summarize experimental data and their SM predictions by using the 
input values in Table [T| We note that all the SM predictions are in good agreement with 
the current experimental data, except the ratio R{D*) which will be discussed in later 
section. For each observable, the relevant parameters for the theory prediction in 2HDM 
type III are enumerated. Apparently, those parameters will be constrained by corresponding 
experimental data. The detailed discussions are presented in the following sections. 

As discussed in the previous sections, the relevant model parameters we are interested 
in aligned 2HDM type III include three mass parameters Mjj±, Ma, and four Yukawa 
couplings Arr, Afcfe, Att, and Act. Here, we choose the light neutral Higgs boson h as the 
observed Higgs boson at the LHC and adopt the alignment limit |^0fH3] . For other choice 
that the heavy neutral Higgs H is observed one, we refer to Ref. [111113]. Direct searches for 
charged Higgs bosons have been performed at LEP [16|, Tevatron [13EH] and LHC [l9l 150] . 
The LEP Collaboration put the lower bound Mh± > 79.3 GeV by assuming B{H^ —)■ r+z/) + 
—)■ cs) = 1 within 2HDM [IH]. The neutral Higgs search at the LEP experiment also 
put lower bound on the neutral Higgs masses such as Mh > 92.8 GeV and Ma > 93.4 GeV 
within GP-conserving MSSM scenario [51]. We adopt those lower limits for heavy Higgs 
masses as reference values even though above results may depend on Yukawa structure and 
msusY scale. Indeed, the lower limits of Higgs masses are irrelevant to our main result. 
With all these considerations, we restrict the parameters of 2HDM type HI in the following 
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observable 

SM 

EXP 

Ref 2HDM parameters 


B{B tv) ■ 10^ 

0.85 ±0.14 

1.14 ±0.22 

[53] A^g, Xifd-) A^it, 

Mh± 

R{D) 

0.297 ±0.017 

0.391 ± 0.041 ± 0.028 

|54| {Xbb), Xtt, Xct, 


R{D*) 

0.252 ±0.003 

0.322 ±0.018 ±0.012 

[Mj (Afefe), Xrr, Xct, 


Amrf[ps“^] 

0.51 ±0.06 

0.510 ±0.003 

[53] (Xbb), Xtt, Xct, 


Ams[ps“^] 

16.93 ± 1.16 

17.757 ±0.021 

[53] {Xbb), ^bs, ^tt, Xct, 

B{B Xs-f) ■ 10^ 

3.36 ±0.23 

3.43 ±0.22 

|S3| Xbb, \t, Xct, Mh± 


B{t-^ eg) 

< 10-^° 

< 1.6 X 10-4 (95% CL) [55] (Xbb), Xu, Xct, (M^^±), 

Mh, Ma 

a{pp tt) 

- 

< 62fb (95% CL) 

|56| Act, Mh, Ma 


Rb 

0.21576 ± 0.00003 0.21629 ± 0.00066 

jSZ] {Xbb), Xtt, Xct, Mh± 


Po 

1 

1.00040 ±0.00024 

|58| Mh± , Mh , Ma 



TABLE II. SM predictions and experimental measurements for the observables used in the numeri¬ 
cal analysis. The last column denotes their dependence on the 2HDM parameters. The parameters 
in the parenthesis imply that they can be safely neglected. 


ranges: 

Mh± e [80,1000]GeV, 

Mh (Ma) e [ 125 (93), 1000] GeV. (8) 

These choices of parameter regions are shown to be reasonable in later section. 

In order to derive an allowed parameter space, we impose the experimental constraints 
in the same way as in Refs. I121E2]: for each point in the theoretical parameter space we 
span the range of the theory prediction for an observable by performing the 2a variations 
of inpnt parameters. If the difference between the central valnes of theory prediction and 
experimental valne is less then the snm of two errors in qnadratnre, then this point is 
regarded as allowed. Since the main theoretical nncertainties are dne to the hadronic inpnt 
parameters, common to both the SM and the 2I1DM, the relative theoretical nncertainty is 
assnmed to be constant at each point in the parameter space. 


IV. TOP QUARK FCNC PROCESSES AT COLLIDERS 

The LHC is often called top-factory since the top pair is copionsly prodnced throngh QCD 
interaction. The LHC Rnn I data already collected millions of top pair events, and even 
mnch more top pair events are expected to be collected in the LHC Rnn H. Undonbtedly, 
the LHC provides us unique chance to explore the top quark FCNC processes which are 
extremely small in the SM. 
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The experimental search for top quark FCNC can be performed either by anomalous 
decays or production of top quarks at hadron colliders with top quark FCNC couplings [5^ 
EU. We note that the searches for f —)■ ch |65l [66] do not provide any constraints on 
2HDM type III in alignment limit since the top quark FCNC couplings with the SM Higgs 
vanish. The anomalous top decays via t —)■ c/uV where H = 7 , Z are explored at the 
Tevatron [671169] and at the LHC [TOllTH] . without finding any significant excess of signal 
events. However, these searches do not provide any meaningful constraints on 2HDM type 
HI since the prediction is much suppressed by loop correction and EW couplings. Contrary 
to top decays, the anomalous single top production has much chance to probe top quark 
FCNC coupling due to the large gluon luminosity in the parton-distribution-function (PDF) 
and the relatively large QCD coupling. The experimental searches for single top events put 
upper bound on B(t —)■ eg) and B{t —)• ug) [711178] . We focus on B(t —)■ eg) by ignoring 
M-quark involved FCNC process since it is extremely suppressed in Cheng-Sher ansatz even 
though u quark PDF is bigger than c quark PDF. 

The same sign top pair production is a tree-level process and therefore promising to 
test NP scenarios which contain top quark FCNC couplings. Notable example is that the 
NP scenario with Z' mediated top quark FCNC coupling [Tsiisn] that explains the anoma¬ 
lous top forward-backward asymmetry observed at the Tevatron [8T1183] is disfavored by 
non-observation of the same sign top pair production at the LHC [HU ES] • The recent ex¬ 
perimental search at ATLAS with integrated luminosity of 20.3 fb~^ at 8 TeV puts the most 
stringent upper limits on a{pp —>■ tt). We interpret the result as an upper limit on ee —)■ tt 
process to constrain A^. 

In what follows, we study the phenomenology oi t ^ eg and ee —)■ tt processes within the 
2HDM type HI to investigate the top quark FCNC coupling. 


A. t ^ eg 

In the SM, t ^ eg decay is extremely suppressed due to GIM mechanism. However, this 
rare top decay can be enhanced in some NP scenarios [HSl EZ]. In general, the form factor 
for the effective teg vertex is dehned by [27] 

tG^T^ ( 9 ) 

IbTT^ \ rrit rrit q^ hJ 

where T“ (a = 1,..., 8) denote SU{3) generators. The form factors A, B, C and V have 
been calculated in various types of 2HDM P |HHl EH]. In the 2HDM type HI, these form 
factors are generated by the penguin diagrams mediated by the neutral Higgses h, H and 
A and charged Higgs H^. Their explicit expressions are given in Appendix]^ With the 
convention Eq. (§. the decay width for t ^ eg is given by m 

^ In Ref. [27], the last two terms of Eq. (j^ are omitted. Although they do not contribute to the width 
r(t —>• eg), they are necessary to satisfy Ward identity. 
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FIG. 1. Branching ratio of t —cgf as a function of the charged Higgs mass. Dashed line: a common 
scalar mass Mh± = Mh = Ma is taken. Shaded region: neutral Higges’ masses Mh and Ma vary 
but constrained by the oblique parameter Ap. 


with Cf = — l)/2iVc. We note that B(t —>■ eg) is proportional to as can be 

seen from Eq. (A2). 

The LHC search for anomalous single top production is performed by ATLAS Collabora¬ 
tion with 14.2 fb~^ at 8TeV [55]. Non-observation of signal put an upper limit on B{t —)■ eg) 
as 

B{t ^ eg) < 1.6 X 10-\ (11) 


In Fig. we show the plot of 2HDM type III prediction for B{t ^ eg) as a function of 
the charged Higgs mass by setting XctXtt — 10 The shaded region is spanned by changing 
neutral Higgses masses under the constraints from Ap. We refer to Ref. |13| for detailed 
analysis of Ap. Even though there can be up to factor O (10^) enhancement comparing to 
the SM expectation for the small Mh±, the current experimental bound is far above the 
theory prediction. Therefore, it would be hard to constrain the top quark FCNC parameter 
space with anomalous single top production measurement at the LHC. 


B. cc —)> tt 

The same sign top pair production at hadron collider requires FCNC coupling with t— 
or M— channel exchange of neutral particle with spin 0 or 1 since the electric charges of 
hnal states are same. Another possibility is s-channel process mediated by a charge 4/3 
new particle. Various NP scenarios that contribute to the same sign top pair production 
are well summarized in Ref. [90] with effective operator formalism. The production rate of 
the same sign top pair at hadron colliders via the contact interactions with different chiral 
conhguration is modeled in Ref. [HI]. Meanwhile, in this work we perform the full theory 
analysis with spin 0 Higgs boson as a mediator since the effective operator formalism may 


^ Our numerical result is consistent with Fig. 3 of Ref. m by setting ^ct = 6t = 1- 










FIG. 2. (a) Total cross section for cc —>■ tt at the LHC 8TeV run in {Mh,Ma) plane. We set 
Xct = 15. The shaded region (green) is allowed parameter space at 95% CL. (b) The allowed 
parameter space in {Mh{= Ma), A^) plane in the case where H and A are degenerated in mass. 


not reproduce well the full theory result if the mediator mass is quite less than 1 TeV. We 
refer to Ref. [^2] for the analysis with another mediators. 

In the 2HDM type III with alignment limit, the same sign top pair production arises at 
tree level via t- or u-channel diagrams with exchange of heavy neutral Higgs bosons, H or 
A. The partonic scattering cross section for qq —)■ tt process is described as 

a(s) = j i) + gA{s, i) + ^intf(s, i)^ , (12) 


where the amplitude square functions Qi are dehned as 


g^{s,t) = KC, 


2^4 
cl 


t — m 


+ 


u — 


+ 


tu — m^s — ml 


^intf (s, i) = 2A^cC 


ct 


(tu + mfs — mf) (tu + (Mjj + M\)(s/2 — ml) + 
(t - Ml)(t - Ml)(u - Ml)(u - Ml) 


(13) 


where cf) = H, A. Then the total cross section for cc —)■ tt is convoluted with parton lumi¬ 
nosity function fcc(x-it‘'F) of sea quark pair cc as follows 


a(cc —)■ tt) 


dxa(xs)fcc(x,tiF ), 


(14) 


where r = Aml/s and fcc(x^g>F) is defined by 


fcc(x, t^F) 



fc/p(y 1 i^f) fc/p(x/y^ Hf) ■ 


(15) 


Here, fc/p(y, A^f) is c-quark PDF and the factorization scale g.F is set to be /ij? = m^. We use 
MSTW2008LO PDF set [93] for the numerical analysis. The gluon and charm quark initial 
state process with extra jet radiation is not considered by assuming that the contribution is 
subleading. 

The experimental searches for the same-sign dileptons and 6-jets at CMS with 19.5 fb~^ [M] 
and at ATLAS with 20.3 fb“^ [56] at 8 TeV can be applied for constraining the same-sign 
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top pair production rate. The non-observation of any significant excess of signal events sets 
the upper bound of the production cross section. The strongest bound comes from ATLAS 
result. ATLAS provides different upper bounds depending on the helicity conhguration of 
effective operators within contact interaction model. We conservatively adopt the largest 
upper bound among the three as follows: 


a{pp —)■ tt) < 62fb (ATLAS 95% CL [56]) . 


( 16 ) 


We re-interpret this result to constrain the cross section a{cc —>■ tt) using the formula 
described above. The constraint is usually strong for small Higgs masses. Since the signal 
rate is proportional to the large values of Act are severely constrained and conversely 


the small value of Act is hardly excluded. Fig. 2(a) shows the prediction of scattering cross 


section by setting Act = 15 in {Mh, Ma) plane and the allowed region with shaded green 
color. As shown, the interference effect is constructive. For the given Act value the region 


MhiMa ^ 400 GeV is excluded. Fig. 2(b) shows the allowed parameter space in {MH,Xct) 
plane for the case where H and A are degenerated in mass. Experimental bound provides 
quite stringent upper limit on Act as 10 ~ 20, depending on the heavy Higgs mass. 


V. FLAVOR PHYSICS - TREE-LEVEL PROCESSES 

Since the top-quark FCNC couplings take part in charged Higgs Yukawa sector, they can 
contribute to the semi-leptonic decay and leptonic decay of B mesons which are tree-level 
processes. In this section we study the two r-involved tree-level processes, B —)■ 
and B ^ TV io constrain top quark FCNC couplings. The former (latter) is involved with 
b —)■ c{u) charged current. Therefore, any NP model which contains such charged current 
with a new charged particle can contribute to these processes [951 - 198] . 

For those processes with b —>■ c{u) charged current, the effective Hamiltonian is described 
by [22] 


BeS — CvLL^VLL + ^S^RL^SRL + ^SLL^SLL^ {Q — '^i c) (17) 

with the effective four-fermion operators 

<^VLL = 

Olu. = (qPRb){fPL,pA 

OLl = (qPLb)(fPRl >,). ( 18 ) 

Within the SM, the vector boson W~ is exchanged, therefore only (Pvll generated with 
tree-level Wilson coefficients 


^<j,SM _ 
'-'VLL — 


AGpVqi, 


( 19 ) 


where Gp denotes the Fermi coupling constant and Vqh are the CKM matrix elements. On 
the other hand, within the 2HDM type HI the scalar charged Higgs boson is exchanged, and 
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therefore C^|ll ^srl generated. The corresponding tree-level Wilson coefficients are 


a 


c, 2HDM 
SLL 


^tb^ct^T 

~Mi: 


^ g, 2HDM 
^SRL 


^gbibbi r 


( 20 ) 


We neglect which is proportional to Xut and extremely suppressed by w-quark mass. 

For B —)■ TV decay, we can dehne a theoretically clean observable by taking the 

ratio with relatively clean measurement B —)■ = e,/i,r) to cancel the hadronic 

uncertainties: 


R(D<-y ). 


( 21 ) 


B{B D(*)iv) ■ 

Note that the CKM matrix element Wfe is also canceled out. Then, the theory uncertainty of 
R{D^*^) are very small, 6(1)%, while the experimental error is quite large, 12(7)% because 
of missing neutrino in r reconstruction. 


With the effective Hamiltonian in Eq. (17), the theoretical prediction of R{D^*^) relative 


to the SM value is described as p^MlnT]. 


R{D) — i?gM(71) j 1 -|- l.SRe 


R{D*) = Rsm{D*) l + 0.12Re 


^SRL + C^SLL 

+ 1.0 

^SRL + ^SLL 

1 

^c,SM 

'-'VLL 

^c, SM 
'-'VLL 



/^c 

'-'SRL 


/-ic 

'-'SLL 


^c, SM 
'-^VLL 


+ 0.05 


r^c 

'-'SRL 


r^c 

'-'SLL 


^c, SM 
'-'VLL 


( 22 ) 


Due to the spin of D* meson, the NP effects on R{D*) are much smaller than the ones on 


R{D) |100i 110211104] . The relevant Wilson coefficients are given in Eqs. (19) and (20). Since 
Cg^L is suppressed by rrib/v in Cheng-Sher ansatz and also by CKM matrix element, its 
contribution is negligibly small. 

The BaBar experimental data for B —)■ have shown somewhat large values com¬ 

paring with the SM expectations for both R{D) and R{D*) where the combined discrepancy 
was 3.4 ct level [251 1105] . It was also discussed that these can not be simultaneously accom¬ 
modated by 2HDM Type II. To explain both discrepancies it was shown that the large top 
quark FCNC coupling which contributes to Cg^LL i'^ E*!- (22) is needed [23l El]. Very 
recently, the Belle collaboration reported the measurements of both R{D) and R{D*) [106] . 
and the LHCb collaboration did for R{D*) [107] . Even though the Belle result is in the 
middle of the SM expectation and the BaBar result, due to the reduced errors, the average 
values are still in 3.9 ct discrepancy [5l| (See Table 2 for comparison). 

The allowed parameter space in (A^-t, A^) with different charged Higgs masses constrained 
by R(D^*'>) is shown in Fig. For any given charged Higgs mass both Xct and Xrr do not 
simultaneously become zero. For small A,-,- value, X^ must be very large. Interestingly, larger 
charged Higgs mass requires larger X^- These feature can be understood as a whole since 


only the product XctXrr/M^^ enters the contributions from 2HDM, as show in Eq. (20). 
Explicitly, the current B —data put the bound 


-0.0030 < XctXrr/M^± < -0.0023, 


(23) 
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FIG. 3. Constraints on and A^r from R{D) and R{D*). Pink-colored region is allowed at 95% 
CL in (Att, Act) plane with different charged Higgs masses. 


which can be seen in Fig. It is noted that A,-,- is associated with the nentral Higgs decay 
H/A —)■ TT. If Xrr is large, the LHC has a good opportnnity to detect nentral Higgs bosons 
in their tanonic decay channels. In the case of small A^r, the conpling Act should be large, 
which may be severely constrained by the same sign top pair production as shown in previous 
section. 


Contrary to i? —)■ decay, B ^ tv decay is a helicity suppressed process and 

more strongly suppressed by CKM factor. Therefore, B —)■ {[ile)v decays are extremely 
rare, 0(10“^) and 0(10“^^) respectively, and not yet measured although B ^ fiv will be 
measured soon at Belle H. Thus, we have no way to cancel the large theory uncertainty of 
hadronic current of H —)■ tv. The uncertainties from the SM prediction and experiment 
for B{B —)■ tv) are very large, 24% and 19% respectively. Due to these large errors, the 
constraint from B{B ^ tv) is not much significant. 


With the effective hamiltonian in Eq. (17), the branching ratio oi B ^ tv reads 


B{B —)■ tv) 


Stt 


2 

-mlmBTBfl 



1 + 


m\ 


^s’rl ~ G, 


SLL 


mbirij 


a 


n,SM 

VLL 


( 24 ) 


where Jb denotes the H-meson decay constant. The relevant Wilson coefficients for 2HDM 
type HI are shown in Eqs. (19) and (20). We note that not only but also ^bs,^bd can 
contribute to Cgj^L within Cheng-Sher ansatz due to the relatively large CKM factors. Even 
^ut can signihcantly contribute to CgLL. Due to the combination of these contributions 
to a single observable B{B —>■ tv), none of these Yukawa couplings get any meaningful 
constraints. 
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VI. FLAVOR PHYSICS - LOOP-LEVEL PROCESSES 


A. Bd Xs'j 

As for the loop-induced process we first consider Bd —)■ Xgj decay. Taking the normaliza¬ 
tion with B{Bd —)■ XcCUf,), the dominant theoretical uncertainties from ml and CKM factor 
are canceled out. The effective Hamiltonian for the Bd —)■ Xs'y decay read |1U81 llU9j 

n,s = -^v;,VaJ 2 CiO., ( 25 ) 

^ i=l 

where the explicit expressions of the tree or penguin operators Oi-g can be found in 
Ref. m- The magnetic penguin operators, Oj and Os, which are characteristic for this 
decay, are dehned as 

^7 = Og = (26) 

where mi, denotes the 6-quark mass in the MS scheme, and e {qs) is the electromagnetic 
(strong) coupling constant. The heavy degrees of freedom from the W~ boson contribu¬ 
tion [IIIHII 9 ] and charged Higgs contribution |12nfll22j are integrated out at mw scale, 
and we obtain the Wilson coefficients = ttivk)- They evolve into /i = scale by 

renormalization group equation and consequently resum the large logarithms in perturba¬ 
tive QCD to all order |12311125] . The higher order correction at fi = mb scale should be 
necessarily done [nsHnn]. 

The compilation of all those calculation for B{Bd —t W 7 ) reached at next-to-next-to- 
leading-order (NNLO) in perturbative QCD |130m32] . (For a recent review, we refer to 
Ref. |133] .i For given NP contributions to the theory prediction for B{Bd —t W^y) at 
NNLO is given by [I33] 

B{Bd Xs'y) X 10^ = (3.36 ± 0.23) - 8.22ReC^P - LOOReCg^P , (27) 

where the hrst number represents the most up-to-date SM prediction. By using current 
experimental data, we obtain 

8.22ReC'^P + l.OOReC'^^ = -0.07 ± 0.32 . (28) 

Therefore, it is natural for to become (P(O.l). 

In the 2HDM type HI, the one-loop contribution to via charged Higgs exchange is 
described by 


^2HDM 

'-^7,8 



(29) 
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where the loop functions -^7 8^^ are given in Ref. [120] and xw = Yukawa 

components and normalized by SM ones are defined as 



It should be emphasized that the A^ term is enhanced by the spin-flip factor mt/rrib and 
becomes comparable to Au. Therefore, it is unique for —)■ that the coupling can 

be significantly constrained. Another interesting feature is that the coefficient Xct of second 
factor in Au is highly suppressed while the one in first term contains CKM-enhanced factor. 
The Xct prefers to be 0(1^ from B Thus, the Xu and Abb must be strongly 

correlated to satisfy Eq. (28). In order to avoid large cancelation between l/SAb^F^g and 
XbbFjg in Eq. (29) that causes fine-tuning, we prefer to take the region where Xbb,Xtt ~ 


0 ( 0 . 1 ). 

To be more specific regarding the fine-tuning argument, we refer to Ref. |134j and re-define 
fine-tuning parameter A for an observable as follows 


max(5(5j) 


( 31 ) 


Here, Q denotes the difference between theory prediction and experimental data and 6Qi 
represents each individual contribution of the theory to the Q. Therefore, small means 
significant fine-tuning. (For example, A = 25 correspond to 4% fine-tuning.) The allowed 
parameter space in {Xu, Abb) plane for given X^ = 10 and Mh+ = 400 GeV is shown in Fig. 
by requiring A“^ > 10%. The gray region causes significant fine-tuning. We note that by 
avoiding significant fine-tuning, not only Xu is constrained but also Abb is highly restricted 
as we expected. 


B. - Bd,s mixing 

The Bq — Bq{q = d, s) mixing occurs via box diagrams by exchanging W boson or 
charged Higgs within 2HDM between Bq and Bq. We note that the tree level diagrams 
can also contribute through h-s-{H/A) vertices within 2HDM type HI. We first study the 
NP contribution from loop processes while the tree-level contribution is discussed in the 
next section. The mass difference Am^ between the two mass eigenstates and B^ is 
related with off-diagonal element of mixing matrix Ml^ such that Am^ = 2|Mj'2l- Since the 
constraints from Bd — Bd mixing appears to be more or less weaker than those from Bg — Bg 
mixing, we only consider latter one in this work. The effective Hamiltonian with AB = 2 
for the Bg — Bg mixing is described by [133 

Y. ( 32 ) 
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FIG. 4. Yellow region is allowed parameter space at 95% CL from 13{Bd —)• Xg^) with fixed Xct = 10 
and M[j± = 400 GeV by requiring > 10%. The gray region causes significant fine-tuning. 


In the SM, only operator can contribute, where 


{FYPls^) ■ (33) 

The corresponding Wilson coefficient is C^^^{mw) = ^Sq{xw) where xw = The 

function S'o(a:) can be found in Ref. [lOSj . Then the Am^ is obtained as 

Am, = 2\{B,\n^^=^\B,)\ = ^^\v:,Vts\^flBBmB^Vhml,S,{xw) • (34) 

Here, = 0.552 is a short-distance QCD contribution. As for the long distance non- 

'^ 1/2 

perturbative quantity Jb^B^^ , we use Lattice QCD result. 

Within the 2HDM, two additional operators are generated by the box diagrams with 
charged Higgs boson exchanged: 


^SRR ^ (6“P^s“) {b^Pns ^), {b^a>^^PRS ^). (35) 


Using the formulae in Ref. |136j . the corresponding Wilson coefficients are obtained as 




xh± + 1 2,xh± log xh± 


_ {xh± - 

_{xh± 

+ 


1)2 {xh± 
-4 -|- X]y 


+ 


1)3 

{xw - ^XH±)\ogXH± 


l){xw - 1) {xh± - IY{xh± - Xw) 
3xw loga^m 


(xw - ^)‘^{xh± - Xw) 




2 2 


/m2 \ 

2 

{xh± + l)logX/b±' 

X^w) 

1 

1 

to 

1 

cm 

1 

Ti 

1 


(36) 


where xh± = The subscript WH or HH represent the exchanged particles in the 

box diagram. We note that = 0 at the matching scale f^w- Contrary to the B^ —)■ Wy, 
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FIG. 5. Allowed parameter space at 95% CL by Anis experimental data for fixed = 500 GeV 
in (a) (Att, Act), (b) (Att, Xct^tt) and (c) (Re[C'y'^'^], Im[C'y'^'^]) planes. Green color (SI) corresponds 
to the solution without significant fine-tuning. Black color (S2) and gray color (S3) represent the 
parameter space with significant fine-tuning, < 10%, where S3 causes large ImM ^2 while S2 
does not. The dashed (dot dashed) line denotes 68% CL (95% CL) bound from The red point 
represents the SM prediction. 


the Ad contribution in has signihcant suppression factor thus its contribution 

is negligible. Although the operators and are generated through operator mixing 
during renormalization group evolution as described in detail in Refs. |1351I137T - I141j at NLO 
QCD, the effects are minor and we do not include them. 

Therefore, only A„ is numerically relevant in Bg — Bg mixing. It contains Xu and Act as 
dehned in Eq. (30) which are constrained by experimental data of Arug given in Table [n} 
The allowed region for the parameter space in (Att, Act) plane as well as (Att, XctXu) plane 
are shown in Fig. |^(a) and (b). We perform more detailed study on the allowed parameter 
space by considering the hne-tuning argument to £t the data. As shown in Eq. (30), there 
are two solutions for = 0 which give the result consistent with experimental data: 


Att ^ ^ (2.14 - 0.04i)Act, 

Vts V m 

V* frrT 

or Att —Act--0.004 Act. (37) 

The parameter space near these two solutions are allowed, but can cause signihcant hne- 
tuning. We represent the allowed parameter without signihcant hne-tuning, or > 10% 
by green color, and for < 10% by black color. 

In the region where the signs of Act and Xu are same, the two 2IIDM contributions 
and destructive with each other. The parameter space that brings the cancelation 

between the two can be another solution to ht the data, but also causes signihcant hne- 
tuning. We represent the parameter space near the solution with signihcant hne-tuning, 
A“^ < 10%, with gray color. For this solution space, the real parts of the two 2IIDM 
contributions are strongly canceled, but sizable imaginary parts still remain as can be seen 
in the Fig. [^c). This sizable imaginary part can cause large time-dependent GP-asymmetry 
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phase in 6 —>■ c decays from the relation 0^“ = arg(Mf 2 ). We show the bonnds at 68% 
and 95% CL in Fig. |^c) with current average value [S3] 

= -0.015 ± 0.035 . (38) 

As shown, the gray region is excluded by at 68% CL, but survives at 95% CL. This 
region will be more signihcantly covered by future experimental data. 

For later convenience, we summarize the features of each parameter regions and their 
color notation with the dehnition of SI, S2 and S3 as follows 

51 : (green color) > 10%, 

52 : (black color) < 10%, ~ 0, 

53 : (gray color) A“^ < 10%, ReC^^j:^ + — 0 , large ImM *2 • (39) 


VII. COMBINED ANALYSIS AND FUTURE PROSPECT 


We hrst combine the constraints from —)• Bg — Bg mixing, and cc —?• tt on the 

couplings Act and Xu- We also include the constraints from EW precision measurements, 
Z ^ bb and Ap. We refer to Ref. [13| for the details of these EW precision measurements. 
We scan the parameter space as described in Eq. ([^. The allowed parameter space is 
obtained by requiring that it accommodates all the experimental data with 95% CL. The 
result is shown in Fig. |^a) for Mh± = 500 GeV. As discussed in previous section we divide 
allowed parameter region into SI, S2 and S3 whose features are portrayed in Eq. (39). 

For the region SI, the requirement A“^ > 10% in Bg — Bg mixing gives the upper bound 
on Act and is slightly stronger than the one from cr(cc —)■ tt) combined with Ap. The 
upper bound on Xu for the region SI is given by Rf,. On the other hand, for the regions 
S2 and S3, the couplings Act and Xu are bounded by a{cc — tt) accompanied with Ap 
and Rb- Therefore, the same sign top pair production plays crucial role to constrain Act 
regardless of hne-tuning. But if we avoid signihcant hne-tuning (for SI), Bg — Bg mixing 
put the signihcant bound. The projection for the exclusion limit at 14TeV with 300 fb~^ is 
estimated by assuming that the statistical error is dominant (See Ref. [14211143] for more 
details about the projection method). The result is outstanding. The upper bound of X^ 
reach 8~15 with 300 fb“^ at 14TeV as shown in Fig. We note that Bd —?• Xg'y does not 
put bound on Xu nor Xu for any parameter sets due to sizable contributions from Xbb term. 

We turn to the B —)■ decays. With hxed Xrr-, B decays also put 

bounds on Mh± and Act. By taking A^r = 40, the allowed parameter space is shown in 
blue-colored region in Fig. (with Mh± = 500 GeV). As shown in Fig. |^b), |Act| has 
different upper limits for each parameter set depending on Mjj±. They lead to lower limits 
on |Arr| as can be seen in Eq. ( [^ and Fig. The allowed parameter spaces in {Mh±, 
|A^t-|) plane are presented in Fig. For hxed Mh±, the lower bounds for S2 and S3 are 
same and slightly diherent from SI. It should be noted that these lower bounds become 
stronger as Mh± increases. Gonversely, the Mh± is upper bounded when A^r is hxed. In 
the case of relatively heavy charged Higgs, the lower bound on Xrr is very strong. With the 
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FIG. 6. Combined constraints from Bg — Bg mixing, cc tt, Z ^ bb and the oblique parameter 
Ap on the 2HDM parameters. The allowed regions are divided into three parts and shown in the 
green (SI), black (S2) and gray regions (S3), (a) Allowed parameter space in {Xu, Xct) plane for 
the fixed Mjj± = 500 GeV. The constraints from cc —)• tt, Ap and Z ^ bb are shown in dashed and 
dot-dashed lines respectively. The projection for cc —?■ tt, Ap at 14TeV with 300 fb~^ data is shown 
by a dotted line. The allowed parameter space by i? —)■ (with A^r = 40) are indicated by 

the blue region, (b) Allowed parameter space in {Mu±, X^t) plane. Note that the upper and lower 
bounds of black region are same with gray region so they are not shown in the plot. 


constrains of cc -X- tt at 14TeV with 300 fb~^ data, the lower bonnd on A^r would become 
twice of current bound as shown in Fig. For Mh± > 500 GeV, the coupling A,-,- should 
be greater than 30, which can signihcantly enhance H/A —)■ tt decays. Therefore, this can 
be constrained by heavy Higgs search with tt final states at the LHC. However the signal 
strength of gg — H/A —)■ tt process strongly depends on heavy Higgses masses and is 
effectively proportional to A^^. Since there are much parameter space near Xu ~ 0 in the set 
SI (green region) as shown in Fig. |^that may avoid the constrains from gg -x- H/A —)■ tt, 
the constraints would be restricted. Perhaps, some part of parameter space, especially small 
Act and large Xu, Xtt region will be excluded. On top of that, for such very large r Yukawa 
coupling, the perturbativity would be threatened. 

We now discuss about the constraints from t —)■ eg. With the above allowed regions SI, S2 
and S3, we make theoretical predictions for B{t —)■ eg). Since the combined constraints put 
upper bounds on both Act and Xu, Therefore, XctXtt is upper bounded in all three parameter 
sets SI, S2 and S3. Note that the set S3 represents also the lower bounds for both Act and 
Xu that comes from Z ^ bb and cc —)• tt as shown Fig. 6(a). The upper bound of H(f —)■ eg) 
for SI, S2 and the allowed region for S3 are presented as a function of Mh± in Fig. 

The current LHC upper limit is much larger than these theory predictions. Thus, it does 
not give any constraints. The projection for the upper limit at 14TeV with 300 fb~^ data is 
also drawn in Fig. in dotted line. As shown, it would be hopeless to see or constrain the 
top quark FCNC couplings from the t ^ eg measurement. 

So far, we have neglected the tree-level contribution to Bg — Bg mixing through the down- 
type FCNC couplings b-s-{H/A) with the Yukawa coupling ^hg. Even though ^i,s is severely 
suppressed in Cheng-Sher ansatz such as ^bs/^bs = 3.6 x 10“^, the tree level contribution with 
0{\) Xgb has no CKM suppression, and is comparable to the loop contribution. By including 
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FIG. 7. Allowed parameter space |At-tI) by the combined constraints from loop-induced 

processes and B —)• decays. See Eq. (39) for the definition of each parameter set. The 

projection for cc —)> tt,Ap at 14TeV with 300 fb“^ data is shown by a dotted line. The lower 
bound of black region is same with gray region so is not shown in the plot. 



FIG. 8. 2HDM prediction on B{t ^ eg) as a function of Mu±. The green and black lines denote 
the upper bounds from the solutions of the combined constraints SI and S2, respectively. The gray 
region corresponds to the solution S3. The dashed line denotes the current upper bound at LHC, 
while the dotted line is for the future sensitivity at 14TeV with 300 fb“^ data. 


the tree level contribution, the allowed parameter space in (A^t, A^) plane is signihcantly 
extended since the large NP contribution from the loop processes can be canceled by the 
tree level contribution. Therefore, including the tree level contribution in Bg — Bg mixing 
always weakens the constraints on Xu, Xct- To understand the effect of A^f, quantitatively, we 
show a plot in Fig. [^for allowed region of Xct Xu with respect to the hxed Xgb value by imposing 
the constraints from cc tt, Z ^ bb and Ap. We see that for Xgb > 0.003 Mh,a large XctXu 
is required to cancel the large tree-level contribution. In fact, for Xgb — O.OOSMha, fhe 
magnitude of tree-level contributions is already comparable to the magnitude of the SM 
contributions. For Xgb < 0.003 Mh,a , the bound on X^Xu is not much changed from the one 
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FIG. 9. Allowed parameter space in {Xgb/M, XctXtt) by the combined constraints from Bg — Bg 
mixing, cc tt, Z ^ bb and the oblique parameter Ap. Mh = Ma = M and Mij± = 500 GeV 
are taken. For Bg — Bg mixing, both the tree-level and loop-level contributions are included. 


M^±=500GeV, Mh=Ma=M 



0.002 0.004 0.006 0.008 


Xsb/M [GeV-‘] 


given in previous section. 


VIII. CONCLUSION 

The general 2HDM as an extension to the SM is a potential NP candidate. To avoid 
severe constraints from down-type quark FCNC, we adopt Cheng-Sher ansatz. This NP 
scenario permits presumably large top quark FCNC coupling Xct, which is the main target 
to be explored in this work with collider phenomenology as well as flavor constraints and 
EW precision measurements. To this end, we consider anomalous single top production 
which can limit —)■ eg) and the same sign top pair production via cc —)■ tt at the LHC 
in association with not only flavor tree-level processes, B —)■ B ^ tv but also flavor 

loop-level processes, B^ —)■ Bg — Bg mixing. 

We hnd that among them the B —)■ D^*Wu, Bg — Bg mixing and cc —)■ tt play important 
role to constrain Xct- Especially, still large value of Xct is preferred by average value of 
measurement with the new data for B —from Belle and LHCb. To bring 
solid understanding of the result, we separate the allowed parameter space into three sets, 
SI, S2 and S3, regarding the hne-tuning to ht the data and the features reflected in the 
observables of Bg — Bg mixing. SI does not suffer from the hne-tuning while S2 and S3 
cause signiheant hne-tuning to ht the data. More specihcally, S3 shows large imaginary part 
of M ^2 while SI and S2 do not. 

For the allowed parameter sets SI, S2 and S3, Xct is severely upper-bounded by either 
cc —)■ tt or Bg — Bg mixing. Therefore, to ht the values, the Yukawa coupling X-rr 

is lower bounded for given charged Higgs mass Mh± and conversely Mh± is upper bounded 
for hxed Xrr- The large Xrr will be constrained hj gg ^ H/A ^ tt, however it strongly 
depends on neutral Higgses masses and Xu- The extended study with heavy Higgs search 
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data at the LHC can be a future work. Since XctKt is small for all the parameter sets and 
the theory prediction is loop-suppressed, the upper limits for —)■ eg) do not provide 


constraints on the remaining parameter space with current experimental data nor in future 
LHC experiment. On the other hand, large Xct is mostly constrained by cc —)■ tt process 
regardless of hne-tuning. cc —)■ tt would play more important role to probe top quark FCNC 
at the LHC 14 TeV Run. 
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Appendix A: Form factors in t ^ eg 

In general 2HDM, the form factors for teg vertex was hrst calculated in Refs. [2711144] . 
Here, we recalculate these form factors and write them in terms of scalar one-loop functions. 
Each form factor in Eq. (|^ is summation of four different contributions from the penguin 
diagrams with A, H and exchanges, e.g. A = Aa + Ah + Ah±- They are calculated as 



VA = -gs^^f^, VH = gs^Hf2, VH± = -gs\Vtt\^^H±ft ■ (Al) 


To compare with Refs. [271 1144] . we neglect the small term in tbH^ vertex of Eq. ([^ 


and show the result in general with complex Yukawa couplings 





(A2) 


The loop functions are dehned as 
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/f = m?(C + Ci^ + Cf), 

ft (A3) 






The scalar one-loop functions are abbreviated as 




m^, 0 , m‘^, m^, 


C'p* = 0,0,0, M^±), 


(A4) 


which are dehned in Refs. |145H147] and can be numerically evaluated by the LoopTools 
package |147] . In the penguin diagrams with charged Higgs we omit the terms propor¬ 
tional to ^bb as in Refs. since these terms are suppressed by nib/v. In addition, we 

have analytically checked that the form factors presented in this paper are in agreement with 
those obtained in Ref. 123 except one minor discrepancy: for the parameter dehned in 
Ref. [27], we obtained + (1 — x)M‘fj But this does not come into play in our 

numerical analysis. 
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